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The  exact  detection  probabilities  and  false-alarm  probabilities  for  continuous 
square-law  detection  of  a  signal  in  Gaussian  noise  are  calculated,  with  no  restric¬ 
tive  assumptions  about  the  size  of  the  observation-time-bandwidth  products.  Two 
cases  are  considered:  they  are  called  the  low-frequency  (LF)  and  narrowband  (NB) 
cases.  In  both  cases,  the  signal  bps  coherent  ana  inconerent  components. 

Q/  _ 

Results  for  the  detection  probability  versus  signal-to-nolse  ratio  (S/N),  with 
false-alarm  probability  as  a  parameter,  are  given  for  both  tlie  LF  and  NB  cases 
for  observation-time-bandwidth  products  equal  to  ^  (n  =  0  (1)  7)  and  for  the  frac¬ 
tion  of  coherent  signal  power  equal  to  0,  .  2,  1,  -For  a  specifiea  detection  proba¬ 
bility  of  .  5  and  false-alarm  probabilities  of  (jW"“  (n  =  1  (1)  8),  the  required  S/N  i.s 
also  plotted  for  the  same  range  of  time-bandwidth  products  and  fractions  of  co¬ 
herent  signal  power  as  above. 
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ABSTRACT 

The  exact  detection  probabilities  and  false-alatm  probabilities  for  square-law  detection  of  a 
signal  in  Gaussian  noise  are  calculated,  with  no  restrictive  assumptions  about  the  size  of  the 
observation-time-bandwidth  products.  Two  cases  ztc  considered:  In  the  first,  called  the  low- 
frequency  (LF)  case,  the  signal  is  a  combination  of  a  dc  component  and  a  Gaussian  process  com¬ 
ponent  with  exponential  correlation.  Here  the  noise  is  a  Gaussian  process  with  identical  exponential 
correlation.  In  the  second  case,  called  the  narrowband  (NB)  case,  the  signal  is  a  combination  of  a 
sine  wave  component  and  a  Gaussian  process  component  with  an  exponentially  damped  cosinu¬ 
soidal  correlation  at  the  same  center  frequency  as  the  sine  wave.  Here  the  noise  is  a  Gaussian 
process  with  identical  exponentially  damped  cosinusoidal  conelation.  The  two  components  of  the 
signal  are  called  the  coherent  and  incoherent  components,  respectively,  in  both  the  LF  and  NB 
cases. 

The  output  of  the  square- law  detector  is  integrated  continuously  over  the  observation  time 
and  compared  with  a  threshold  for  statements  about  signal  presence.  It  is  not  assumed  that  the 
square-law  detector  output  is  sampled  infrequently  enough  to  yield  independent  samples. 

Results  for  the  detection  pi obability  versus  signal-to-noise  ratio(S/N),  with  false-alarm  proba¬ 
bility  as  a'  parameter,  are  given  for  both  the  LF  and  NB  cases  ior  observation-time-bandwidth 
products  equal  to  2”  (n  =  0(1)7)  and  f  -  the  fraction  of  coherent  signal  power  equal  to  0,  .2,  1. 
For  a  specified  detection  probability  of  .5  and  false-alatm  probabilities  of  10*"  (n  =  1(1)8),  the 
required  S/N  is  also  plotted  for  the  same  range  of  time-bandwidth  products  and  fractions  of  co¬ 
herent  signal  power  as  above.  Comparisons  with  a  Gaussian  approximation  for  large  time-band¬ 
width  products  are  made;  it  is  shown  that  the  Gaussian  assumption  is  optimistic  in  predicting 
performance  of  the  scuare-law  detector,  the  exact  amount  depending  on  the  time-bandwidth 
product. 
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OPERATING- CHARACTERISTICS  fOR  CONTINUOUS 
SQUARE-LAW  DETECTION  IN 
GAUSSIAN  NOISE 


1.  INTRODUCTION 


Ofteiii  signals  in  noise  are  detected  by  filtering  a  received  process  to  the 
frequency  band  of  interest,  squaring  the  filtered  output,  continuously  integrat¬ 
ing  over  an  observation  .intery^  where  the  signal  is  anticipated,  and  comparing 
this  decision  variable  with  a  threshold  or  with  the  past  history  of  this  variable. 
In  fact,  under  certain  conditions,  the  optimum  detector  for  stochastic  signals 
takes  this  form  (for'  example.  Ref.  1,  Chapter  XI).  This  form  of  processor 
could  be  employed  whether  the  received  signal  is  deterministic  (known  or 
unknown),  stochastic,  or  a  combination,  of  both,  although  it  is  not  necessarily 
optimum  in  all  these  cases.  This  processor  is,  in  fact,  frequently  adopted  for 
passive  detection  of  targets  of  xmknowh  statistics. 

I  ' 

Two  important  par^eters  that  characterize  the  performance  of  this  proc¬ 
essor  are  the  bandwidth  of  the  stdchastic  process  at  the  input  to  the  squarer 
I  and  tbe  .observation  (or  integration)  time.  Practical  cases  occur  where  this 
time-bandwidth  product  ranges  from  the  order  of  unity  to  several  thousand.  In 
the  latter  case  of  very  large  time-bandwidth  p.roducts,  the  decision  variable 
can  often  be  accurately  approximated  by  a  Gaussian  random  variable  (RV). 
However,  when  the  time-bandwidth  product  is  of  the  order  of  unity,  it  can  not 
be  approximated  as  Gaussian,  especially  when  very  low  Mse-alarm  probabil¬ 
ities  or  Iiigh-detection  probabilities  are  of  interebt.  The  question  then  arises 
as  to  the  exact  performance  of  the  processor.  Calculation  of  the  performance 
is  complicated  by  the  difficulty  in  analyzing  the  continuous  integration  of  a  non- 
Gaussian  process  since  a  "sum”  of  dependent  random  variJibles  is  formed.  Also, 
approximation  of  the  performance  by  using  a  sum  of  independent  variables  is 
tenuous,  because  the  exact  value  of  the  "effective"  number  of  independent 
samples  to  use  is  not  ea.sily  as  certained.  Thus,  there  is  a  need  for  (1)  the  exact 
evaluation  of  the  receiver  operating  characteristics  (ROC)  of  the  square-law 
detector  followed  by  an  integrator,  both  as  a  benchmark  against  v/hich  other 
processors  a-^d  approximations  c^  be  compared,  and  (2)  the  determination  of 
when  the  Gaussian  assumption  can  be  validly  employed,  especially  for  medium 
time-bandwidth  products  andivery  small  false-alarm  probabilities. 
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Two  cases  of  signal  and  noise  spectra  at  the  input  to  the  squarer  are 
considered:  In  the  low-frequency  (LF)  case,  the  signal  is  a  combination 
of  a  dc  component  and  a  stationary,  zero-mean,  Gaussian  process  component 
with  exponential  correlation,  and  the  noise  is  a  Gaussian  process  (of  different 
power  level)  with  identical  exponential  correlation.  In  the  narrowband  (NB)  ease, 
ti?.e  signal  is  a  combination  of  a  sine-wave  component  and  a  Gaussian  process 
component  with  an  exponentially  damped  cosinusoidal  correlation,  and  the  noise 
is  a  Gaussian  process  with  identical  correlation,  except  for  power  level.  The 
two  components  of  the  signal  are  called  the  coherent  and  incoherent  componen'  s, 
respectively,  in  both  the  LF  and  NB  cases,  for  ease  of  discussion. 

Actually,  the  general  derivations  of  processor  performance  capability  pre¬ 
sented  here  are  not  restricted  to  the  cases  above.  They  actually  allow  for  an 
arbitrary  deterministic  coherent  signal  component  (which  may  be  known  or  im- 
known)  and  an  arbitrary  (but  identical)  correlation  function  for  the  incoherent 
signal  component  and  the  noise.  However,  numerical  evaluation  of  performance 
capability  requires  evaluatin.Y  jh<:!  resolvent  (or  reciprocal)  kernel  and  the 
Fredholm  determinant  of  thr  correlation  function.  Hence,  numerical  results 
have  been  confined  hers  to  the  cases  given  in  the  paragraph  above. 

Several  practical  applications  fall  under  the  framework  above.  For  example, 
if  the  incoherent  component  of  the  received  signal  in  the  NB  case  is  absent,  we 
have  detection  of  a  sine  wave;  this  is  encoimtered  in  active  or  passive  detection 
of  a  target  or  in  binary  on-off  communication  with  another  party.  Or,  if  the 
ceberent  component  is  absent  in  the  LF  or  NB  case,  we  have  detection  of  a 
random-signal  process  in  noise;  this  situation  can  hold  true  in  long-range  trans¬ 
mission  through  a  fluctuating  medium  or  in  reflection  of  an  acoustic  signal  off 
the  time-varying  ocean  surface.  It  is  often  employed  in  passive  detection  of  a 
noisy  source.  Finally,  if  both  signal  components  are  present,  relevant  examples 
in  the  NB  case  are  reflections  of  an  acoustic  signal  from  a  surface  with  low 
roughness,  or  communication  via  a  medium  supporting  both  a  direct  (nonfluc¬ 
tuating)  path  and  fluctuating  paths.  Extensions  to  multitone  applications  are  also 
possible  and  ;u’e  discxissed  in  Section  4. 

A  great  deal  of  work  on  the  statistics  of  thetime-averagepowerofa  Gaussian 
process  has  been  presented  in  the  past.  Rice^  gave  the  first  four  semi-invariants 
(or  moments)  for  an  arbitrary  correlation  fimction  of  the  Gaussian  process, 
and  approximated  the  probability  density  function  (PDF)  for  very  small  and  large 
observation  tiitnes.  Emerson’  gave  an  infinite  product  for  the  characteristic 
function  (CF)  of  a  filtered  squared- Gaussian  process  and,  via  the  first  three 
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semi-invariants,  approximated  the  PDF  through  Gram-Charlier  and  Laguerre 
expansions.  Siegert'*  ga%'e  the  CF  in  terms  of  the  solution  of  an  integral  equa¬ 
tion  obtained  by  Kac  and  Siegert.^  In  the  former,  the  CF  was  la  terms  of  an 
infinite  product,  although  it  was  shown  that  the  resolvent  kernel  offered  an  al¬ 
ternative  approach.  No  examples  were  presented.  Darling  and  Siegert®  and 
Siegert’  showed  that  the  CF  of  a  more  general  functional  could  be  obtained 
from  the  solution  of  two  integral  equations,  and  they  obtained  the  CF  in  closed 
form  for  a  finite-tirae  exponential  averagei; for  a  Gaussian  process  with  an  expo¬ 
nential  correlation.  Slepian®  obtained  a  series  expansion  for  the  PDF  and  cumu¬ 
lative  distribution  function  (CDF),  where  each  term  was  given  by  an  integral. 
No  extreme  probabilities  were  evaluated  and  the  accuracy  was  about  1  percent. 
Jacobs^  considered  a  communications  problem  with  random  signals  and  made 
a  finite  approximation  to  the  CP  by  increasing  the  number  of  terms  imtil  the 
CDF  essentially  stabilized.  Schwartz  obtained  a  closed  form  CF  for  correla¬ 
tion  functions  whose  corresponding  homogeneous  Fredholm  integral  equation 
reduces  to  an  equivalent  linear  differential  equation.  Four  examples  were  given, 
and  both  LF  and  NB  cases  were  considered.  Steenson  and  Stirling  ”  conducted 
a  simulation  and  gave  results  on  the  PDF  and  faise-alarm  rate  and  showed  that 
the  gamma  distribution  accurately  fitted  the  data.  Some  related  results  on  the 
distribution  of  quadratic  forms  are  given  in  Deferences  12,  13,  and  14. 

None  of  the  previous  results  are  of  sufficient  accuracy  to  allow  calculating 
the  cumulative  distributions  very  near  zero  or  unity  --  the  regions  of  interest 
for  very  small  faise-alarm  and  high-detection  probabilities.  The  complicated 
results  for  the  CF,  even  when  in  closed  form,  indicate  the  need  for  quick  and 
accurate  CDF  calculations  that  do  not  require  moment  evaluations  or  series 
expansions.  Here  we  make  no  attempt  to  find  analytic  expressions  for  the  PDF 
or  CDF,  but  proceed  directly  nmnerically  from  the  CF  of  the  decision  variable 
to  the  CDF  via  a  Fast  Fourier  Transform  (FFT),'^'**'  The  assumption  that 
the  incoherent  component  of  the  received  signal  has  the  same  correlation  (and, 
hence,  spectrum)  as  the  noise  in  both  the  LF  and  NB  cases,  except  for  power 
levels,  enables  evaluating  the  CF  in  closed  form. 

If  the  integrator  following  the  squarer  were  replaced  by  a  sampler  and  sum¬ 
mer,  calculation  of  the  CF  in  closed  form  can  be  accomplished  fairly  simply 
once  the  eigenvalues  of  the  matrix  of  sampled  correlations  are  evaluated.  This 
technique  allows  very  general  ?nd  different  signal  and  noise  spectra  and  will  be 
the  subject  of  a  future  report.  For  now,  however,  we  restrict  ourselves  to  con¬ 
tinuous  integration  without  sampling  the  squarer  output. 
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2.  PROBLEM  DEFINITION  , 

Two  different  cases  will  be  considered  in  tMs  report,  n^ely,.  theXF  and 
N3  processors.  Block  diagrams  of  each  processor  v.lll  be  given  and  the  perti¬ 
nent  characteristics  of  the  signal  and  noise  described.  Then  the  expected  per¬ 
formance  of  the  two  cases  will  be  briefly  compared  for  the  case  of  large  time- 
bandwidth  products. 

2. 1  LOW-FREQUENOY  CASE 

A  block  diagram  of  the  LF  processor  is  given  in  Fig,  1,  where  y(t)  is  a 
LF  waveform  composed  of  three  components: 

y(t)  =c(t)  +  s(t) -f  n(t)  , 

In  (1), 

a^  c(t)  is  a  LF  deterministic  waveform  (known  or  unknown)  and  is  called 
the  coherent  component  of  the  received  signal. 

b.  s{t)  is  a  stationary,  zero-mean,  LF  Gaussian  process  \vith  correlation 

*^8  e2cp(-2W^  Irl) 


and  Is  called  the  incoherent  component  of  the  received  signal. 

The  noise  n(t)  has  the  same  statistical  description  as  s(t),  but  the  two  processes 
are  independent.  Its  correlation  is 

Rn(n=c2exp(-2Wjr|)  .  (3) 


The  (double-sided)  spectra  of  s(t)  and  n(t)  are  given  by 

G  (0=^2  - - - 

/nl^  ’ 


(4) 


where  W*  can  be  interpreted  as  the  effective  (or  statistical)  baidwidth*"  of  the 
positive-frequency  components  of  the  LF  random  processes: 


w  = 

e 


(0) 


2/  drR2^„(r) 


(5) 
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The  idontical  spectral  shape  of  the  iiiicdierent  component  s(t}  and  the  noise 
a(t)  is  realized  in  practice  by  filtering  fairly  broadband  received  signal  and 
noise  processes  throng  a  low-pass  RC  filter.  It  cw  be  interpreted  physically 
as  overresolving  the  received  signal  vdih  a  very^low^paSs  filter. 

The  decision  variable  i  in  Fig,  1  is  given  by 

£  =4  /  +  sW  +  n(t)1^  (6) 

T  T  J, 

and  is  to  be  compared  with  a  threshold  for  detecting  signal  presence.  (No  as¬ 
sumption  about  the  size  of  the  time-bandwidth  product  TWj  is  made. )  This 
is  broadband-energy  detection  and  contains  three  subcases  of  importance: 

a.  c(t)  =  0.  Signal  presence  indicated  by  an  increased  power  level  in  the 
received  random  process. 

b.  s(t)  =  0.  Signal  presence  indicated  by  a  nonzero,  time- varying  mean  in 
the  received  random  process, 

c.  Both  c(t)  and  s(t)  nonzero.  Signal  presence  indicated  by  adding  a  non¬ 
zero,  time-varying,  mean-random  process  to  the  received  random  process. 

2. 2  NARROWBAND  CASE 

A  block  diagram  of  the  ND  processor  is  given  in  Fig,  2.  In  the  figure,  y(t) 
isan  NB  waveform  composed  of  three  components,  as  in  (1).  The  signal  c(t)  is 
an  NB  waveform  centered  at  frequency  f^ ,  with  a  deterministic  complex  enve¬ 
lope  (known  or  unknown),  and  is  called  the  coherent  component  of  the  received 
signal.  Both  s(t)  and  n(t)  are  independent,  stationary,  zero-mean,  NB  Gaussian 
processes  with  correlations 

(It  is  assumed  that  foT  »  1,  where  T  is  the  integration  time  in  Fig.  2. )  By 
definition,  a(t)  is  called  the  incoherent  component  of  the  received  signal. 


(8) 


The  (double-sided)  spectra  of  s(t)  and  n(t)  are  given  by 


s,n' 


A 

1 

1 

{2-/  . 
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where  can  be  interpreted  as  the  effective  bsudwidlli  of  the  pdsitive-frequeut^ 
components  of  the  NB  random  processes: 


W  —  ^  ^  ^  ^  StO 

2/drR2^{r) 


The  identical  spectral-’shape  of  the  incoherent  compontiiit  s{t)  and  the  noise 
n(t)  are  realized  in  practice  by  filtering  fairly  broadband  received  signal  and 
noise  processes  through  an  RLC  filter  that  is  KB.  It  can  be  interpreted  physically 
as  overresolving  the  received  signal  with  a  very  NB  filtcir. 

The  decision  variable  Jl  is  given  by 

t  ^  /  di  !£(t)  »•  s{t)  +  n  (t)|  2  ,  (10) 

where  x(t)  denotes  the  complex  envelope  of  x(t)  (see  Appendix  A).  The  decision 
variable  JL  is  compared  with  a  threshold  in  order  to  detect  sij^nal  presence. 
(No  assumption  about  the  size  of  the  time-bandwidth  product  TW,  is  made. ) 
This  is  NB  energy  detection  and  contains  three  subcases: 

a.  c(t)  =  0.  Signal  presence  indicated  by  an  increased  power  level  in  the 
received  rajcdom  process. 

b.  s(t)  =  0.  Signal  presence  indicated  by  an  NB  waveform  in  the  received 
random  process. 

c.  Both  e(t)  and  s(t)  nonzero.  Signal  presence  Indicated  by  adding  an  NB 
waveform  and  random  process  to  the  received  random  process. 

2.  3  PERFORMANCE  COMPARISON  BETWEEN  LOW- 
FREQUENCY  AND  NARROWBAND  CASES 

The  effective  bandwidths  of  the  LF  and  NB  processes,  which  have  been  de¬ 
fined  in  (5)  and  (9),  respectively,  are  measures  of  the  spectral  width  of  the 
positive-frequency  components  of  the  respective  processes.  A  bandlimited  low- 
pass  process  confined  to  frequencies  (-W,  W)  must  be  sampled  at  least  eveiy 
(2W)~‘  seconds  in  order  not  to  lose  any  information.  Thus,  in  along  time  interval 
T,  there  would  be  (2W)T  samples  upon  which  to  base  a  decision  as  to  signal 
presence  or  absence.  On  the  other  hand,  a  bandlimited  NB  process  confined  to 
frequencies  (f^  -  W/2,  f^  +  W/2)  and  the  corresponding  negative  frequency  band 


TS  4233 


must  have  its  Li-phase  and  quadrature  components  sampled  at  least  every  W“^ 
seconds  Tin  order  not.to  lose  any  information.  Then  in  a  time  T,  there  would  be 
2(W)T  samples  upon  which  to  base  a  decision  as  to  signal  presence.  Therefore, 
both  bandiimited  processes,  LF  and  NB,  yie?d  the  same  number  of  information 
samples  when  W  is  interpreted  as  the  spectral  width  of  the  positive-frequency 
components.  As  a  result,  in  our  investigation  of  nonbandlimited  processes,  we 
expect  that  when  the  effective  positive-frequency  bandwidths  are  equal,  that  is, 
Wc  =  We ,,  the  performance  of  the  LF  and  NB  processors  should  be  comparable 
for  large  TW^ .  This  will  be-  borne  out  by  the  numerical  examples  in  the  next 
section. 

For  small  time-bandwidth  products,  however,  the  performance  of 
the  LF  and  NB  processors  differ.  For  example,  as  T — -0  in  (6)  and  (10), 
■^Lp“*'-yLP  [y2  (to)  +  y^  (tj]  (see  Appendix  A).  Thus,  in  the 

limit  of  small  time-bandwidth  products,  LF  processing  contains  one  degree  of 
freedom,  v/hereas  NB  pi’ocessing  contains  two  degrees  of  freedom.  The  per¬ 
formance  of  the  processor  will  be  better  than  that  of  the  LF  processor  for 
small  time-bandwidth  products. 
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The  CFs  for  the  LF  and  NB  cases  are  defiv(3d  in  Appendix  B,  the  general 
results  being  (B-22)  and  {B-46).  lliese  results  are  specialized  to  the  cases  of 
interest  here,  resulting  in  (B~31)  and  (B-51),  respectively.  For  the  special 
case  of  T  — ►O,  (6)  and  (10)  reduce  to  a  single  sanaple  of  the  input  processes, 
and  derivation  of  the  detection  probabilities  and  false-alarna  probabilities  is 
possible  in  closed  form;  these  derivations  are  presented  in  Appendix  C. 

The  CDF  P  of  a  nonnegative  RV  from  its  CF  f  can  be  evaluated  numerically 
via  use  of  any  of  the  equivalent  forms**^  : 

F(X)  .i  -if  expUCTI 

4f  (ll, 

=  i-|  ;  ^f.(f)cos(a),  X>o  . 

In  (11),  f,  and  f;  are  therea),  and  imaginary  parts,  respectively,  of  CF  f.  The 
error  in  terminating  the  integrals  in  (11)  at  finite  limits  for  computer  evalua¬ 
tion  is  considered  in  Appendix  D. 

3. 1  LOW-FREQUENCY  CASE 


Suppose  i  in  (6)  is  compared  with  a  threshold  X  for  purposes  of  signal 
detection.  We  now  define  a  normalized  threshold, 


and  the  time-bandwidth  product. 


(12) 


/3  =  TW^  . 


(13) 


Values  of  the  required  normalized  threshold  An  for  different  values  of  0  and 
specified  false-alarm  probability  Pp  are  given  in  Table  1. 


Before  presenting  results  for  the  detection  probability,  it  is  convenient  to 
define  several  power  signal-to-noise  ratios  (S/Ns).  We  define  the  received  co¬ 
herent  signal  power 


T  T 


(14) 
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Table  1 


THRESHOLDS  FOR  LF  PROCESSOR 


1 

0 

1 

1 

2 

4 

8 

16 

32 

r —  - 

64 

128 

2. 706 

2. 065 

1. 841 

1. 634 

1.461 

1.329 

1.233 

1.  264 

1.115 

6.635 

4. 302 . 

3.362 

2.600 

2. 065 

1.710 

1.478 

1. 324 

1.223 

10. 828 

8. 706 

4. 971 

3. 564 

2. 630 

2. 047 

1.684 

1.455 

1.307 

15. 137 

9. 180 

6. 631 

4. 540 

3. 182 

2. 364 

1.  872 

1.571 

1.381 

10-’ 

19. 511 

11.  693 

8. 322 

5. 526 

3,  728 

2. 670 

2. 050 

1.  678 

1.448 

23.928 

14.280 

10.031 

6.523 

4.272 

2. 970 

2. 221 

1. 779 

1. 511 

28. 374 

16.785 

11.753 

7.527 

4.815 

3.265 

2.387 

1.  877 

1.  570 

10~® 

32.  841 

19. 352 

13.483 

8. 537 

5.389 

3.  558 

2.549 

1.  627 

Tna  coherent  and  incoherent  S/JNfs  are  defined,  respectively,  as 


n 


and 


The  total  received  S/N  is  defined  as 


R 


Rc^R. 


(15) 


(16) 


(17) 


We  also  define  the  fraction  f  of  total  received  signal  power  that  is  coherent  as 


(19) 


using  (15),  (16),  and  (17).  Then  we  can  express 

R^  =  fR, 

R^=(l-f)R  , 

and,  thereby,  utilize  the  total  received  S/N  R  and  the  fraction  f  of  coherent 
signal  power  as  the  fundamental  parameters  characterizing  the  performance 
of  the  processor. 

In  the  series  of  Figs.  3  throu^  11,  to  follow,  plots  of  the  detection  prob¬ 
ability  Pp  versus  rV^  are  giv-an,  with  Pp  as  a  parameter.  Each  figure 
corresponds  to  a  different  value  of  0,  specifically  0  -  0,  1,  2,  4,  8,  16,  32,  64, 
and  128.  Each  figure  is  composed  of  three  parts:  Part  A  is  for  f  =  1,  part  B  is 
for  f  =  .  2,  and  part  C  is  for  f  ~  0. 

The  equation 

dli  =  10  log  ,0  R  =  20  log  R  (20) 

canbe  employed  to  convert  the  abscissato  decibels  in  Pigs.  3  through  11.  (R  is 
a  power  S/N,  whereas  is  a  loaeasure  of  a  voltage  S/N. )  Equation  (20)  is  the 
total  input  S/N  expressed  in  decibels. 

The  difference  in  performance  for  different  fractions  f  of  coherent  signal 
power  is  very  marked  for  small  0,  Thus,  Fig.  3,  for  0=0,  illustrates  a  sig¬ 
nificant  drop  in  Pp  for  large  values  of  R’/^  when  f  is  decreased  from  unity.  In 
fact,  this  degradation  in  performance  occurs  even  when  f  is  decreased  from  1 
to  .  9  (this  curve  is  not  shown).  As  8  increases,  the  difference  in  performance  is 
lessened  until,  for  large  0,  the  opposite  behavior  occurs  for  large  rV^  (for 
example,  P^  «  .  9999).  That  is,  the  required  value  of  R’/^  for  a  specified  Pp 
near  unity  is  less  for  f  =  0  than  for  f  =  1  for  very  large  /?;  see,  for  example. 
Fig.  li. 

On  the  other  hand,  for  small  values  of  Pp ,  the  smaller  values  of  f  yield 
the  better  performance  for  small  0,  whereas  the  larger  values  of  f  give  bet¬ 
ter  performance  for  large  0.  There  is,  however,  rather  slight  differences  in 
performance,  that  is,  required  S/N  R,  in  the  case  of  small  Pp. 

The  reason  for  tliis  cross-over  of  the  curves  is  explained  by  the  fact  that 
we  are  considering  a  fixed  receiver  and  varying  the  transmitted  signal  composi¬ 
tion.  The  fixed  receiver  is  more  nearly  optimum  for  some  transmitted  signals 
than  others,  and  it  is  impossible  to  order  the  performance  of  the  various  com¬ 
binations  of  coherent  and  incoherent  signal  power. 
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The  situation  is  made  more  striking  when  the  perfonnance  of  the  optimum 
receiver  for  completely  coherent  reception,  f=l,  is  considered.  For  a  dc 
signal  waveform,  LF  processing,  and  exponential  noise  correlation,  we  can 
use  the  first  equation  on  p.  114  of  Reference  1  to  evaluate  Helstrom’s  power 
S/N  d^  for  the  optimum  receiver  as 

d^  =  R^(l4.Tff^)  =  R{l+TW^)  =  R{l+|8)  .  (21) 

The  false-alarm  and  detection  prdsabilities  for  this  case  ai'e  given  by  [Ref.  1, 
p.  106,  Eq.  (3.32)] 

P_=<I'(-x), 

^  _  (22) 

Pj,  =  <I)(d-x)  =  +j8  \/r-s)  , 

where*  x  is  selected  for  prescribed  Pp .  (For  a  specified  performance  (Pp » 
Pp)  of  this  optimum  processor  for  f=l,  required  values  of  d  can  be  read 
from  Ref.  1,  p.  88,  Fig.  IV.  1. )  When  (22)  is  superposed  on  part  A  of  Figs.  3 
through  11,  it  lies  everywhere  above  the  square-law  detection  curves,  as  is 
expected.  Furthermore,  the  slope  of  the  straight  lines  ((22);  is  greater  than 
the  slopes  of  the  square-law  detection  curves.  However,  when  (22)  is  super¬ 
posed  on  parts  B  and  C  of  Figs.  3  through  11,  the  equation  lies  below  the  square- 
law  detection  curves  for  small  R  and  above  it  for  the  larger  values  of  R.  The 
square-law  detector  can  give  better  performance  in  some  regions  of  R  (for 
f  <  1)  than  the  optimum  coherent  receiver  (for  f  =  1)  for  the  flat  transmitted 
signal,  because  the  transmitted  signal  (the  flat  waveform)  is  not  the  optimiun 
one  to  transmit  in  this  noise  spectrum.  Rather,  a  deterministic  signal  that 
oscillates  arbitrarily  fast  in  the  interval  T  and,  thereby,  occupies  a  region  of 
low-noise-density  does  arbitrarily  good;  there  is  no  optimiun  signal  unle  s  we 
impose  some  bandwidth  constraint.  This  case  has  not  been  considered  here. 

Two  analytical  results  along  a  related  line  are  possible  for  /3  =  0.  From 
(C-7),  it  is  found  that  the  detection  probability  for  fixed  Pp  and  fixed  R(>0) 
always  decreases  as  f  increases  from  zero  (for  A„>0)  regardless  of  the 
values  of  the  false-alarm  probability  and  S/N.  (An  <  0  and  R  <  0  are  uninteresting 
situations.)  Of  course,  for  large  A„,  that  is,  low  Pp,  the  rate  of  decrease 
is  not  very  large  and  may  not  be  noticeable  from  the  plots.  However,  for  large 
An  andR,  for  example,  R  =  An  -1,  (C-7)  gives 


*<^  is  the  CDF  for  a  zero-mean  unit  variance  Gaussian  RV;  see  (C-3). 
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=  _(2;r)-*/2  exp(-l/2)  [  1 - 
0  \ 


(23) 


which  is  significant. 


For  f  =  1,  we  find  the  same  result  for  .S  =  0;  see  (C-S).  The  slope  can.  be¬ 
come  very  negative  for  large  Ap. ;  in  fact. 


/N./ 

f  =  1 

n 


A 


2(2-e)'-'2 


A  -♦CO 


(24) 


Thus,  for  /3=  0,  performance  of  the  square-law  detector  is  degraded  as  f  in¬ 
creases  toward  1  for  certain  S/N  R  and  normalized  thresholds  A  „ . 

For  large  time-bandwidth  products  /5 ,  it  is  anticipated  that  the  decision 
variable  i  in  (6)  could  be  approximated  by  a  Gaussian  RV.  In  Appendix  E  are 
presented  tie  equations  for  the  detection  and  falSv'i-alarm  probabilities  xmder 
this  Gaussian  assumption;  see  (E-3)  and  (E-4).  In  Fig.  12,  these  equations  are 
plotted  for  /3  =  128,  where  f  =  1,  .  2,  and  0,  respectively.  These  approximate 
curves  are  generally  above  the  exact  curves  given  in  Fig.  11  for  the  same  con¬ 
ditions;  that  is,  the  Gaussian  assumption  gives  optimistic  performance  pre¬ 
dictions  compared  to  the  true  value.  However,  the  deviation  is  not  great  at 
j3=  128. 

For  a  specified  performance  (Pp,  Pp),  the  required  input  S/N  R  decreases 
to  zero  as  T  —  ~  .  In  fact,  from  (E-4),  we  find 


Pp~‘l>(x/^‘(i-A^))  as  )S  cc 


(25) 


where  normalized  threshold  A„  is  selected  to  keep  the  argument  fixed;  that  is, 

\/F(A„-l)  =Th(Pp)  (26) 

is  kept  fixed  at  a  threshold  value  that  depends  on  the  prescribed  Pp.  Then  (E-3) 
yields 

Pd~  ^(\^R  -  Th  (Pp))  as  /3  oo  .  ^27) 

Thus,  the  required  S/N  R  decays  as  T"‘/^  as  T  —  «  ,  rather  than  according 
to  T~*;  this  is  attributable  to  the  square-law  detector  characteristic.  Notice 
that  this  rate  of  decay  of  required  R  holds  regardless  of  the  value  of  f;  that 
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is,  it  is  independent  of  the  fraction  of  each  type  of  signal  component.  Althou^ 
tsi-3  rule  was  deduced  from  the  Gaussian  assumption,  it  holds  in  the  general 
case  because  the  decision  variable  JL  approaches  Gaussian  as  j8  —  «  . 

Figures  4  through  11  may  he  used  to  obtain  the  exact  values  of  the  input  S/N  R 
that  are  required  to  attain-a  specified  performance.  In  parficndar,  for  Pp  =  10*“  , 
n  =  1(1)8,  and  = .  5,  the  required  value  of  R  (in  decibels)  is  plotted  versus 
0  in  Fig.  13.  The  three  parts  of  the  figure  correspond  to  £=  1,  .2,  and.O, 
respectively.  The  rate  of  decay  of  R  with  increasing  /?  is  approximately  /3"V^, 
as  discussed  above.  In  decibels,  the  decay  is  -5  log  ^  . 


The  Gaussian  assumption,  Fig.  12,  gives  values  of  S/K  R  that  are  optimistic 
by  amounts  that  depend  on  the  quality  of  performance  desired.  For  example,  the 
Gaussian  approximation  for  /J  =  128,  Pp  =  10*^  ,  and  Pjj  -  .  5  indicates  a 
required  S/N  of  -3. 05  dB  for  all  *  tyhereas  the  actual  required  S/N  is  -1. 98 
dB,  a  difierence  of  1. 07  dB.  In  fact^  the  -3. 05-dB  S/N  would  yield  a  P^ 
greater  than  10“^  at  /S=  128.  (For  /S=  ^4,  the  difference  in  decibels  increases 
to  almost  1. 5  dB. )  However,  for  Pp  =  .  1,  instead  of  10"®.,  the  decibel  dif¬ 
ference  between  the  exact  answer  and  the  Gaussian  approximation  is  only  .  28  dB 
at  d  =  128.  These  comparisons  have  been  made  here  only  at  the  Pp  =  .  5  level, 
and  conclusions  at  other  levels  can  be  very  different.  Figures  3  through  11  fur¬ 
nish  information  for  required  S/N  at  other  performance  levels  of  interest. 


A  note  of  caution  is  in  order  at  this  point:  From  (C-1),  (C-5),  and  (C-6), 
for  0  =  0,  if  f  =  the  decision  variable  j2  is  given  by 


£  ^(p^z  +  n)^ 


(28) 


and 


Pp  =  $(R  _  A V2)  +  (i,{_R  1/2  _  ai/2)  ^ 

Pp  . 

r  n 


(29, 


^‘This  independence  of  f  for  required  S/N  is  true  only  for.Pp  =  .5;  see  {E-3). 
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Now  for  a  fixed  Pp,  A„  is  fixed,  and  Pp  versus  plots  as  a  straight 
line  on  normal-pi*obability  paper  (for  example,  see  Kgs.  3  tfarqugh  11)  for  all 
values  of  An-  However,  the  RV  j2  in  (28)  is  definitely  not  Gaussian.  Thus', 
a  strai^t-iine  plot  on  noirmal-probability  paper  can  ^not  be  interpreted  as 
corresponding  to  a  Gaussian  RV.  "  *  i  .  » 

* 

3.2  IiSARROWBAND  CASE  ,  .  ' 

Most  of  the  comments  in  the  previous  subsection  are  applicable  here  and 
will  not  be  repeated,  for  the  sake  of  brevity.  Al6'^>the  noffcation  fs  similar,  except 
that  , 

I  . 

j8  =  w^T  ‘  (30) 

Values  of  the  required  normalized  threshold  for  different  values  of  ^  and 
specified  Pp  are  given  in  Table  2.  As  anticipated  in  Section  2.13,  the  thresh¬ 
olds  fo?  large  time-bandwidth  product  $  are  virtually  identical  to  those  for.  the  . 
LF  case  (Table  1).  However,  the  thresholds  for  small  $  are  quite  different. 

I 

I 

Table  2 

THRESHOLDS  FOR  NB  PROCESSOR  ' 


0 

1 

2 

4 

* 

8 

i 

32 

i 

1128 

10"^ 

2. 303 

1. 979 

1.799 

1. 615 

1.454 

1.326 

1. 232 

1. 16^ 

1. 115 

10-2 

4. 605 

3.680 

3. 123 

2. 528 

2.  044 

1. 703 

1.475 

1.  324 

I 

1. 222 

10"^ 

6.908 

5.381 

4.447 

3.423 

2.  593 

2.  036 

1.681 

1.454 

1 

1.307 

lO-'* 

9. 210 

7.  083 

5. 770 

4. 316 

3. 129 

2. 349 

1.  868 

1. 569 

1.381 

10-^ 

11. 513 

8.784 

7. 093 

5. 209 

3. 658 

•  2.  652 

2.  045 

1. 676 

1.448 

1 

o 

13. 816 

10.485 

8.416 

6. 101 

4.184 

2.948 

2. 215 

1.  778 

l'  510 

10-7 

16. 118 

12.186 

9. 739 

6.  994 

.4.709 

3.240 

2.380 

1;  874 

1.569 

10-8 

18.421 

13.  887 

11. 063 

7.  886 

5. 233 

3.  528 

2. 541, 

1.968 

1. 626 
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'  ■  s  : 

'  a;  Figs.  14  throu^  22,  plote  of  yer^,  ares  prese’otM  for  the 

,  KB  processor,  wi&  Pp  as  a  p^ameten,For.  yS.f=  128i.ae-^ve5  j^e  vuia^y 
ideMcal  to  those  for  fke/LB  c^e;.  ii[us,  the  G^siah  apprpMTnation  could  be 
^pHed.here..  ecpial  v^dity  as.  for-  theyLF  cas^,.  AIs6,  :diis-  coi±pb6'rat,es 
did  anticipated  .equality  of  performance  for  the  definittons  of  .  effective  b^dwidth 
adopted  earlier.  The  curves  ediihit  behavior  similar  to  that  for  the.LF  case, 
particularly  the  cross-overs  for  different  values  of  f ; 

*  I  ?• 

The  required  input  S/N  R  to  attain  Pp  =-10~",  n  =  1(1)8,  and  Pp  =  .  5  is 
5  plotted  in  Mg.  23  for  the  NB  case.  The  r^uifed  S/N  is  less  than  the  one,ih'the  LF 
case  for  small  but  is  virtually  identical  to  the  liF  case  tor  ’  large  ^  . 
Required  Values  of  R  for  other  detection  probabilities  are  available  from  Mgs. 
14.1hrou^  22..  ’  ■  . 
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4.  DISCUSSION 


In  the  problem  considered  here,  no  assumptions  about  the  size  o.' the  time- 
bandwidth  product  are  necessary,  and  the  decision  variable  need  not  be  ap¬ 
proximately  Gaussian.  Small  false-alarm  probabilities  can  be  easily  and  very 
accxirately  investigated  numerically.  The  complexity  of  the  closed  form  CF, 
even  for  the  relatively  simple  case  considered  here,  precludes  any  analytic 
solution  for  the  CDF;  nevertheless,  very  accurate  answers  for  the  CDF  can  be 
obtained.  This  numerical  approach  is  recommended  for  other  problems  of  this 
type. 


For  each  doubling  of  the  time-bandwidth  product,  it  has  been  shown  that 
the  required  input  S/N  R  decreases  by  1. 5  dB  (for  large  products)  in  order  to 
maintain  a  prescribed  performance  (Pp,  Pp).  This  result  is  similar  to  that  cited 
for  an  envelope  detector  [Ref.  J.8,  p.  262,  and  Exercises  8.  8  and  8. 9]  and  is 
a  result  of  the  squaring  operation  in  the  receiver.  The  improved  performance 
for  f  =  0  versus  f  =  1  at  low  S/N  has  also  been  observed  for  other  processors 
operating  in  fading  media  [Ref.  18,  pp.  208-9  and  pp.  214-5]  . 

Use  of  a  Gaussian  approximation  on  the  decision  variable  for  large  time- 
bandwidth  products  is  in  error  by  approximately  1  dB  when  the  time-bandwidth 
product  is  100.  For  larger  products,  the  Gaussian  approximation  yields  pro¬ 
gressively  better  estimates,  at  least  for  detection  probabilities  that  are  not 
extreme  and  for  false-alarm  probabilities  that  exceed  10  This  applies  to 
both  the  LF  and  NB  processors  under  the  conditions  investigated  here. 


Some  extensions  of  the  problem  studied  here  are  possible.  They  center 
around  the  ability  to  obtain  the  CFof  the  decision  variable  either  in  a  closedform 
or  in  a  reasonable  amount  of  computer  time.  For  example,  suppose  the  decision 
variable  jg  is  given  as  a  weighted  siun, 


N 

t  ^  2  w  £  , 

k=l  '' 


(31) 


where  the  kth  branch  output,  ,  has  signal  s,,(t)  and  noise  n|,(t)  independent 
of  all  other  branches.  (This  occurs,  for  example,  if  a  bank  of  filters  is  excited 
by  broadband  processes,  and  the  filter  in  the  kth  branch  is  disjoint  with  that  in 
the  jth  branch  for  j  k. )  The  statistical  situation  obtaining  Ln  the  kth  branch 
can  be  any  one  of  the  six  cases  listed  under  (6)  and  (10)  and  can  be  dissimilar 
in  different  branches.  Furthermore,  the  time  of  observation,  the  effective  band¬ 
width,  center  frequency,  and  S/N  can  be  different  in  each  branch.  The  CF 
corresponding  to  the  RV  i  in  (31)  is  given  by 
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! 

N 

ii  =^n^  '  (32} 

where  £,,(5)  is  the  CF  of  the  RV  in  the  kth  branch. 

Examples  that  fit  into  the  framework  of  (31)  include  the  following; 

a.  Multitone  active  detection  or  communicadoa  in  a  time*  varying  random 
mediuni  that  gives  rise  to  the  scatter  component  of  the  received  signal.  (The 
multiple  tones  could  be  transmitted  to  combat  fadj.ng  or  jamming  or  both. ) 

b.  Multitone  passive  detection,  where  the  scatter  component  of  the  received 
signal  may  be  due  to  either  medium  pei’turbations  or  source  instabilities. 

c.  Single-tone  passive  detection  by  comparison  with  a  threshold  established 
from  neighboring  filter  outputs. 

A  compilation  of  performance-capability  characteristics  has  been  presented 
for  continuous  processing  of  both  LF  and  NB  processes.  Howevei*,  they  are 
limited  to  an  identical  exponential  correlation  of  the  random  received  signal 
and  noise  processes;  nonidentical  spectra,  for  example,  of  the  signal  and  noise 
processes  are  disallowed.  The  only  other  continuous  processor  that  yields 
closed-form  expressions  lor  the  CF  of  the  decision  variable  is  an  exponential 
integrator  [Ref.  7,  p.  4l]i,  rather  than  the  finite-time  average  considered 
here.  However,  it,  too,  requires  an  exponential  correlation  function.  The  only 
reasonable  way  to  generalize  to  arbitrary  (unequal)  signal  and  noise  spectra  is 
to  consider  sampled  processors.  Often  this  is  the  practical  situation  anyway. 

The  decision  variable  for  the  square-law  detector  can  thenbe  transformed  into 
a  sum  of  squares  of  independent  Gaussian  variables,  and  the  CF  evaluated  as  a 
finite  product.  If  the  size  of  this  product  is  not  excessive,  direct  numerical 
calculation  of  the  CF  and,  hence,  the  CDF  by  the  methods  in  Refs.  15  and  18 
is  possible.  This  approach  will  be  the  subject  of  future  investigation  on  more 
general  processors  and  spectra. 

In  passing,  it  is  noticed  that  the  detection  problem  posed  here  is  theoretically 
a  singular  one  for  the  purely  incoherent  signal;  that  is,  detection  canbe  accom¬ 
plished  with  zero  probability  of  error  (see,  for  example.  Ref.  1,  pp.  304-306). 

However,  the  optimum  method  of  processing  is  impractical  —  it  would  I’equLre 
precise  knowledge  and  use  of  all  the  eigenvalues  of  the  correlation  kernel; 
equivalently,  second-derivatives  of  the  observed  process  are  required.  The 
physical  impossibility  of  accomplishing  this  in  a  practical  case  forces  us  to 
adopt  a  suboptimum  procedure  similar  to  that  considered  here. 
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Fig.  1.  Low-Frequency  Processor 


TTB14233 


SQUARE -LAW  ENVELOPE  DETECTOR 


Fig.  2.  Narrowband  Processor 
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Appendix  A 

COMPLEX  ENVELOPE  KELATIONS 


Ir  this  appendix,  we  collect  the  complex  envelope  felatiohs  that  are  re¬ 
quired  in  analyzing  the  narrov/band  (NB)  processor.  Consider  a  stationary  zero- 
mean  process  x(t)  passed  through  a  time-invariant  linear  filter,  as  depicted 
in  Fig.  *  -1.  Here  U  (f)  is  the  tmit  step  function: 


U(n  = 


f<oj 


(A-1) 


x(t) 


2U(f) 


x^(t) 


Fig.  A-1.  Single- Sided  Filter 


Thus,  the  output  x^{t)  contains  only  positive-frequency  components.  If  the 
spectrum  of  x(t)  is  Gx(f)i  then,  from  Fig.  A-1,  the  spectrum  of  X4.(t)  is 


Also,  since 


G^^{f)=4U(f)G^(f)  . 


U(f)  =  1  +  sgn(f)  =  1  +  i  [-i  sgn(f)3  , 


we  can  express 

x+(t)  =  x(t)  +  ix^(t)  , 

where  x„(t)  is  the  Hilbert  transform  of  x(t): 

x„(t)  ■=—  /  ar - 

n -oo  t-r 


(A-2) 

(A-3) 

(A-4) 

(A-5) 


Let  the  "center"  frequency  of  (f)  be  f  ^  and  let  the  complex  envelope 
x(t)  of  x(t)  be  defined  as 


x_(t)  =  x^(t)  exp(-i  2;Tfpt) 

Then  the  spectrum  of  the  complex  envelope  is 


(A-6) 
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G  Jf)  =G,  (f+£  )  =4U{f+f  )  G  (f+f  )  , 

'  x+  o  o  X  ,0  >  (A-7) 

and  the  mean  magidtude-squared  v^ue  of  the  complex  ravelope  is 

1^^=  f  df  G  (0  =  4  /  df  G  Yf)  -  2al  , 

-  o  (A-8) 

where  <tI  is  the  power  of  process  x(t). 

Let  the  complex  envelope,  which  is  low-frequency  (LF),  be  represented  in 
terms  of  its  real  and  imaginary  components: 

x(t)  =x^(t) +i  x^(t)  .  (A-9) 

There  follows,  using  (A-4),  (A-6),  and  (A-9), 

x(t)  =x^(t)  cos(27rf^t) -x^(t)  sin(2n'f^t)  ,  (A-10) 

Thus,  x<.(t)  and  Xs(t)  are  the  in-phase  and  quadrature  components  of  x(t). 

Now  if  the  process  x(t)  is  filtered  by  a  linear  filter  with  transfer  function 
H(f),  yielding  (complex)  output  y(t),  it  can  be  shown  that 

~y(t)y(t-r)  =  /  df  exp(i2:rfr)  G  ff)  H(f)  H(-£)  .  (A.-11) 


Applying  this  result  to  Fig.  A-1,  we  obtain 

x+(t)  x^(t-r)  =  0  for  all  r  ,  (A- 12) 

since  U(f)  andU(-f)  are  disjoint.  Hence,  from  (A-6), 

x(t)  x(t -r)  =  0  for  all  r  .  (A-13) 

Therefore,  usii^  (A-9)  and  (A-13), 

Xc(t)  x^(r-r)  =  Re  lx(r)l  Im  ix(r-r)| 

=  ~  [x(r)  +  X*  (t)]  ~(x(r-r)  -  x*  (r-  r)] 

2  -  -  j2  - 

(A-14) 

=  i.[_R  (r)  HR*(r)l  =  -i.lmiR  (r)i 
14  i  i  2  - 

=  -— lm</  df  exp(i2rrfr)  G  (o[  =  -  —  /  df  sin{2?7f!')  G  (f) 
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where  we  have  defined 


Rj^(r)  (r— r)  i  ,  (A-15) 


and  utilized  the  fact  that  G,(f)  is  real  (see  (A-7)).  Now  from  (A-14),  it  may  be 
seen  that  if  and  only  if  the  LF  spectrum  Gx^(f)  is  even  about  f  =  0,  then  the 
processes  Xc(t)  and  Xs(t)  are  uncorrelated  for  allrelative  time  (delays.  In  terms 
of  the  original  spectrum  0,(1),  processes  Xj(t)  and  Xsft)  are  uncorrelated  for 
all  time  delays  if  and  only  if  G,(f)  is  even  about  fo  (in  (0, 2fo)).  In  the  special 
case  where  x(t)  is  a  Gaussian  process  (not  assumed  up  to  ncv/),  Xj(t)  and 
Xs(t)  are  independent  processes  [Ref.  19,  p.  4l]  ,  since  Xc(t)  and  Xs(t)  are 
Gaussian  processes  (being  linear  combinations  of  3c(t)). 

Lastly  (eliminating  the  Gaussian  assumption  again),  by  ,a  procedure  analr- 
ogous  to  (A-14),  I 

x.(t)  X  (t-r)  =  X  (t)  X  (t-r)  =~Re{R  (r)!  =~  /  df  cos (2n^fr)  G  (f)  ..(A-16) 

SS  2-2-00  _  i 

As  a  special  case,  the  powers  are,  usii^  (A-8),  '  ' 

^)='^=Y  fdfG^(f)=(72  .  (A-17) 
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.  i  App^dix.B  ‘ 

'  DERIVATION  OF  CHARACTERISTIC  FUNCTION 
FOR  LOW-FREQUENCY  AND.-NARROWBANID  CASES 

i 

s  *  ,  ■ 

This  appendix  consists  of  jfcwo  parts:  The  first  part  derives  the  characteristic 
fimction  (dF)  of  the  decision  vari^fe  Jl  for  the  Iqw-frequency  .(LF)  case,  as 
given  by  (6)*;  the  secoM  parti  does  the  same  for  the  narrowb^d  (NB)  decision 
variable  i,  ;as  given  by  (10).  *  An  example  for  each  case  is  presented. 


LOW-FREQUENCY  CASE 


From  (6),  * 


i  = —  /  dtrc(t)  + s(t)  +  n(t)]^  . 
T  T 


The  powers  in  the  zOiO-meah  Gaussian  processes  s(t)  and  n(t)  are 
<r  2 ,  respectively.  Let  us  now  define 


i  r(t)=- 


c(t) 


(c'?  ■+  O'?) 


1/2 


(B-1) 
cr2  and 


Mt)  =■ 


s(t)  +  n(t) 


(cr?  +  <T^) 


1/2 


(B-2) 


Then 

!t^(a2+a^)p  ,  (B-3) 

!  I 

where  I 

ps  —  /  dtCr(t)  +  x(t)]^  .  (B-4)  , 

'  ,  T  T  I 

Let  the  correlation  of  x(t)  be  R(r),  which  is  general  here,  except  that  R(0)  =  1, 
as  seen  fi‘om  (B-2);  s(t)  and  n,(t)  are  assumed  independent  of  each  other.  At 
this  point,  the  deterministic  waveform  c(t)  is  also  general. 

*This  equation  is  in  the  main  text  of  this  report. 
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The  CF  iji  of  the  random  variable  is  related  to  the  CF  fp  of  the  RV  P 

according  to 


using  (B-3). 

Now  we  define  for  the  kernel  (correlation)  R,  the  eigenvalues  |  Xj  ^  ,  and 
eigenfunctions  : 

A.  /  du  R(t-u)  ^j(u)  =  4>^t),  t  tT  (B”6) 

T 

Next,  we  expand  x(t)  in  the  set  of  eigenfunctions  according  to 


x(t)  =  S  x.^.(t),  ttT  ,  (B-7) 

where 

X,  =  / dt  x(t)  ,  (B-S) 

since  the  eigenfunctions  are  orthonormal  [Ref.  1,  pp.  98-103].  TheRVs  |xj| 
have  the  properties 


Xj  =  0  , 

_ ;  1  ai-9) 

R(t-u)  ^6.(t)  =  —S 

T  )•< 

using  the  fact  that  x(t)  has  zero-mean  and  employing  (B-8)  and  (B-6).  Since 
X,  in  (B-8)  is  a  linear  operation  on  the  Gaussian  process  x(t),  then  jxj}  aie 
independent  Gaussian  variables. 

When  we  employ  (B-7)  in  (B-4),  we  obtain 

p  =  ~/  dt  r^(t)  +  —  Sr.x.  X?  ,  (B-10) 

T  ■'t  T  i  >  '  T  i  ’ 

where 

r.  =  /j.  dt  t  (t)  ^fi)  .  (B-11) 
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,(^)  =  exp(i^p)=exp^iiy'dtr2(t)^n|e*P  ‘  ‘  ^  [ 


.  (B-12) 


using  (B-10)  and  the  independence  of  the  |  x  j  |  .  The  jth  average  in  (B-12)  is, 
letting  (Tj  =  l/Xj  and  using  (B-9), 


/•  dx 

/  \  / 

/^/27a. 

exp  1  — - 1  exp  1 ! 

\  2a?/  \ 

M  A 

1 - r.  X  +  1 — X 

T  ’  T 


') 


-1/2  /  1,2  jt2/T-2 

exp  (  - 


2rf  en^  \ 


(B-13) 


4-  -i2^T 


) 


)■ 


where  the  square  root  is  principal  value.  Equation  (B-12)  canthenbe  expressed 
as 

_  1  / 

u'2 


— ? 

■p2  i 


r? 

) 


Aj-i2^/T 


•(B-14) 


However,  employing  (B-11),  we  get 

.2 


2L( - ll  uC  du  r(t)  r(u)  2L/ - 

i  A  -i2^T  JJj  i  x.-i 

*  ) 


^j(t)  0,(u) 

2^T 


(B-i5) 


dt  du  r(t)  r(u)  P  (t,  u;  i  2  ^T)  , 


where  F  is  the  resolvent  (or  reciprocal)  kernel  of  R  [Ref.  20,  p.  141,  Eq,  64 
and  p.  134,  Eq.  52] .  Substituting  (B-15)  into  (B-14)  and  employing  (B-2)  yields 

-1/2  . 


= 


2^ 


s  n 


exp 


•  ^  1 


s  n 


cHt) 


(B-iG) 


du  c(t)  c(u)  r  (t,  u;  =  2  ^/T) 
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Now,  if  the  coherent  signal  component  is  absent,  that  is,  c(t)  =  0,  this  CF 
must  reduce  to 


(B-17) 


where  D  is  the  Fredholm  determinant  [see  Ref.  10,  p.  18]  (notice  that  the 
eigenvalues  |Aj|  here  are  equal  to  |l/Xj}  in  Ref. '10).  The  sqitare  root  in 
(B-17)  is  not  principal  value,  but  must  be  determined’  by  tracing  D(i  2  { /T) 
from  ^  =  0,  where  D(0)  =  1,  in  a  continuous  fashion.  This  result  follows  by 
inspecting  the  product  of  principal  value  square  roots  in(B-16).  Utilizing  (B-17) 
in  (3-16),  we  obtain  the  CF  of  p  as 


D(i2^/T) 

-1/2  r 

exp 

- 

2^' 


a 


2 

s 


+ 


dt  du  c(t)  c(u)  r  (t,  u;  i2  ^/T) 


(B-18) 


This  CF  depends  on  the  Fredholm  determinant,  on  the  resolvent  kernel  of  cor¬ 
relation  R,  and  on  the  waveform  c(t).  If  the  Fredholm  determinant  and  the 
double  integral  on  the  resolvent  kernel  can  be  expressed  in  closed  form,  the 
CF  of  p  is  available  in  closed  form. 


It  is  convenient  to  normalize  waveform  c(t)  as  follows:  Denote  the  average 
power  in  received  waveform  c(t)  over  the  interval  T  as 

~  4  >  (B-19) 

and  define  the  normalized  waveform 


(B-20) 


The  average  power  in  Cj(t)  over  T  is  imity.  Also  we  define  two  S/Ns: 


R 


c 


(B-21) 


These  are  the  ratios  of  coherent  signal  power  to  noise  power,  and  incoherent 
signal  power  to  noise  power,  respectively,  at  the  input  to  the  LF  processor 
in  Fig.  1.  Equation  (B-18)  can  then  be  expressed  as 
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fp(^  =[D(i2^/T)]**''^  exp 

"  irt  r(t,u;i2  ^/T)j  . 


(B-22) 


Ex^ple:  Exponential  Correlation;  Constant  Waveform. 

We  assume  that  s  (t)  and  n(t)  are  independent  and  have  the  same  exponential 
correlation,  as  given  in  (2). Then,  using  (B-2), 


R(r)  =e}:p(-2W  1  r|  ) 


The  coherent  waveform  is  taken  to  be  dc: 


c  j(t)  =  1 ,  t  t  T  . 


Then,  from  Ref.  10,  Eq,  (37),t  we  have 


|e(t)  =0 


=  exp  (/3)  - - sinh  +  cosh  (2/3(^^//S)) 


where 


/3  s  w^T, 

^u)  s  (1  -i2u)^^^  . 


(B-23) 


(B-24) 


{B-25) 


{B-26) 


The  square  root  for  <t>  is  principal  value;  that  for  fp  is  not,  but  must  be  traced 
continuously  from  1  =  0,  where  fp(0)  =  1. 

The  resolvent  kernel  is  [Ref.  10,  p.  21,  Eqs.  (29)  and  (30)] 


where 


r(t,  u;i2^/Tj  - - [cos(<«l)  cos((uu)  -  5in(cijt)  sin(<uu) 

cjT  y 

-  y  sin(<uit -u|)]  ,  t, u,  sT, 


cjT  =  i2/30(^/i8)  s  12/395 


(B-27) 


(B-28) 


*  This  equation  is  in  the  main  text  of  this  report, 
t  The  factor  4  in  Ref.  10,  Eq,  (37),  rhould  be  a  2. 
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and 


1  +cosh(2/3^(^/iS))  +  ^(^/(S)  sinh(2/S^(^//3)) 

y  =  -i - : - 

sinh(2/3^(^//3))  +  [  1  +  cosh{2j89i(^//S))] 

cosh(|3^)  +  ^  sinh(^<^) 

=  -i -  (B-29) 

sinh(/3^)  cosh(j8<^) 

I  +  <f>  tanh(j8^) 

^  +  tanh 


There  follows,  after  some  manipulation, 


/Tdtdu  r  (t,u;i2^/T)  = - ^ - - 

^  ''“'T  /3(l-i2f/|8) 


1- 


sinh(2/3 


)3 <^[  1  +  cosh(2/3^)  +  <^sinh(2/3^)]  ) 


(B-30) 


Finally,  the  CF,  (B-22),  becomes 


^  r  1  -  i  ^//3  1 

(^  =  exp(/3)  -  sinh(2i8^)  +cosh(2/3^)J 


1/2 


exp 


Rc 

i - - • 

1  +R. 


^  8  2^  sinh(2/3<^) 

-  <1 -i - 

l-i2^/^  (  |g2^  1 +  cosh(2/39i)  +  0sinh(2/3<;6) 

where  (p  is  given  by  (B-28)  and  (B-26)  as 

9i-(l-i2e//3)*^^  . 

(As  0-^0,  there  follows  from  (B-31) 


fp(^)  -» (1  -  exp 


R. 


1  +R^  1  -12^ 


(B-31) 


(B-32) 

(B-33) 
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which  may  be  readily  derived  as  the  CF  of  p  in  (B-4)  for  one  sample  .f 
[r(t)  +  x(t)]  ^  .  The  CDF  is  derived  in  Appendix  C. ) 

The  CDF  can  be  found  from  (B-31)  by  employing  the  methods  in  R(  1.;.  15 
and  16.  The  square  root  in  (B-31)  is  not  principal  value,  but  must  be  traced 
continuously  from  ?  =  0,  where  fp(0)  =  1. 


NARROWBAND  CASE 


From  (10),  * 


=  f  l£.< 

2T  Jr 


.(t)  +  5(t)  +  n(t)  ] 


(B-34) 


The  powers  in  zero-mean  Gai;issian  processes  s(t)  and  n(t)  are  a\  and  , 
respectively.  We  now  define 

c.(t) 

i(t) 


/  2  2\1/^ 
K  *  <^n> 


(B-35) 


2l(0  = 


s.(t)  +  n(t) 

/  2  2\  *  /  2 


As  a  result, 


where 


+  a2)  p  , 


=  -i-  C  dt  I r(t)  +  XI 


2T^t 


,(t)l 


(B-36) 


(B-37) 


A  A  A 

The  correlation  of  x(t)  is  2R(r),  v/here  R(t)  is  general,  except  that  R(0)  =  1 
(see  Appendix  A). 


The  characteristic  function  (CF)  of  the  random  variable  (RV)i  is 

•  (B-38) 


’^This  equation  is  in  the  main  text  of  this  repoit. 
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In  order  to  evaluate  the  CF  f  p ,  we  expand  r(t)  and  x(t)  in  terms  of  their  real 
and  imaginary  parts: 

l(t)  =  r  (t)  +  ir  (t)  , 

'  '  (B-39) 

i(t)  =  Xg(t)  +  ix^(0  . 


Then 


where 


p  =— (Pe  +  P<) 


tJt 

=  ^f  dtl 

TJt 


[r.(t)  +  x^(t)]  , 


(B-40) 


(B-41) 


[r.(t)+xrt)l  . 


Now  if  and  only  if  the  spectrum  G(f)  corresponding  to  R(t)  is  even  about  f  =  0 
(that  is,  Gs^n(f)  even  about  fo  in(0,2fo)),  then  Xc(t)  and  Xs(t)  are  independ¬ 
ent  random  precedes  [Ref.  19,  p.  4l]  ,  each  with  correlation  R(t).  Then, 
using  (B-40),  the  low-frequency  (LP)  results  of  (B-16),  and  the  symmetry  prop¬ 
erty  r(t,u;X)=  r(u,t;  X)  of  the  resolvent  kernel  of  R  that  is  symmetric, 
we  obtain  the  CF  of  p  as 


f  (^  =  f 


4^^  1  rr 

- II  dt  du  c(t)  c*(u)  r  (t,u;i^/T) 

s  n  J 


The  average  power  in  c(t)  over  T  is,  for  T  »  1, 

=  -L  r  dt  l£(t)  1^  . 

T  Jt  2T  Jt 

We  now  define  a  normalized  waveform. 


(B-43) 


c,(t)  =(2P^)'‘^\c(t)  . 


(B-44) 
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for  then 


T-X  * 


(B-45) 


The  CF  in  (B-42)  can  now  be  expressed  as 
=[D(i^/T)] 

R 


■  r.  ^  , 

exp  1  -  ^ 

L  l  + 

"■c  2  1  rr  •  .  1 

- ^  -  11  *  <1“  £i(t)  C  (u)  r  (t,u;i^/T) 

1  +  R^  J 


(B-46) 


following  the  reasoning  developed  in  (B-16),  (E-l?),  and  (B-18),  and  using  (B-21,. 
Notice  that  this  CF  depends  on  the  complex  envelope  ‘  jCi{t),  and  not  just  on  the 
physical  envelope  |  c  ^  (t)|  of  the  coherent  signal  component  of  the  narrowband 
(NB)  deterministic  waveform.  However,  any  constant  imkno’.vn  phase  in  this 
NB  deterministic  waveform  is  irrelevant  because  it  cancels  out  of  the  expres¬ 
sion  £i(t)£*(u).  That  is,  if  complex  envelope  Cj(t)  is  expressed  as 


Cj(t)  =  E(t)  exp[i^(t) +i0]  , 


(B-47) 


where  E(t)  is  the  envelope  modulation,  <l)(t)  is  tb.e  phase  modulation,  and  d 
is  an  xmknown  phase  independent  o:  ime  t,  then 


(B-48) 


£l(t)  cJCu)  =  E(t)  E(u)  exp  [i<^(t)  -  i^(u)] 

is  independent  of  6, 

Example:  Exponential  Correlation;  Constant  Complex  Envelope. 

Here  we  assume  the  correlation  R(r)  of  complex  envelope  x(t)  to  be 
exponential: 

ft(r)^exp(-^Jr|);/3.^^T  . 

The  correlation  of  the  corresponding  NB  waveform  is  given  by  a  form  similar 
to  (7).  *  We  also  take. for  the  complex  envelope  of  the  coherent  signal  component 


£ , (t)  =  1  , 


(B-50) 


*This  equation  is  in  the  main  text  of  this  report. 


95 


TR  4233 


V/bich  corresponds  to  a  pure  sine  wave  for  c(t).  When  an  approach  analogous 
to  that  given  above  for  the  LF  case  is  employed,  the  CF  (B-46)  becomes* 


f(^=[f  ^/2, /3  -  ^/2)r 


"LP 


1  —i^/0 

=  exp(/§)  - 7 - sbih  (Pd>)  +  cosh(j3  <p)  j  exp 

L  9  J 


R- 


1  +  R. 


•  (B-51) 


1-4 


sinh  (j8 


l-i2^/$  {  1  +  ccsh{/3^)  +  $  sinh  (/3^ 


si] 


where 


(B-52) 


and  we  havfe  utilized  (B-31).  (As  ^  0,  {B-51)  yields 

exp 


(B-53) 

U+R,  1-i^  J 

which  may  be  readily  derived  as  the  CF  of  the  RV  p  in  (B-40)  for  one  sample. 
The  cumulative  distribution  function  (CDF)  is  derived  in  Appendix  C. ) 


The  CDF  can  be  fojind  from  (B-51)  by  employing  the  methods  in  Refs.  15 
and  16.  There  is  only  one  square  root  involved,  that  for  ^  in  (B-52),  and  it  is 
a  principal  value. 


*Tlie  leading  factors  in  (D-51)  agree  with  Ref.  10,  Eq.  (44),  when  we  note  that  his  is 
T/2=^/2  here. 
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DETECTION  PROBABILmES  FOR  ONE  SAMPLE 


As  the  observation  time  T  in  (6)*  and  (10)*  approaches  zero,  the  decision 
variables  A  take  special  forms  for  which  the  distributions  can  be  evaluated  in 
closed  form.  Also,  the  signal  and  noise  spectra  are  not  restricted  in  any  way. 
Here  we  utilize  some  of  the  notation  developed  in  Appendix  B. 


LOW-FREQUENCY  PROCESSmO 
From  (6)*  and  (B-19),  for  T  —  0, 

£  =(P'^  +s+n)^  s  (P^  +  .  (C-1) 

Then 

Prob  >X)  =  ProbdP^  +xl>X) 


since  x  is  Gaussian  vnth  variance  <^1+^1  .  Here  is  the  cumulative  dis¬ 

tribution  function  (CDF)  for  a  Gaussian  random  variable  (RV), 

<t(x)  S  /*  dt(2ff)‘^  exp(-tV2)  ,  (C-3) 

•Ok' 

and  the  normalized  threshold  A„  is  defined  in  (12).  * 

We  define  the  fraction  f  of  the  total  received  signal  power  that  is  coherent 
as 


*This  equation  is  in  the  main  text  of  this  report. 
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»c  •  K 


Pc  +  R.+R.  ^ 


(C-4) 


where  R  is  the  total  S/N  at  the  inpiut  to  Fig.  1.  Then  (C-2)  yields 

>1/2/  I 


p  -d - 1  +$ - ^ — 

\(l+R(l-f))^^/  \(l +R(l -£))*''■ 


(C-5) 


for  the  detection  probability.  The  false-alarm  probability  follows  by  setting 
R  =  0:  •  ‘  ' 

P_  =2$(-A*^)  (C-6) 

F  n  , 

•  I 

Notice  that  the  signal  and  noise  spectra  do  not  appear  in'  (C-5)  dr  (C-6). 

From  (C-5)  there  follows 

**  £  *4  *  I  /*>  «. 

(C-7) 


r  A  '1  A‘/2.R 

- ^  =  -(2ir)'’*''^  exp  I - r;  "-rl  — - 

d{  L  2(l+R)J  (i+R) 

f  =0 


J/2 


Also, 


dPr 


(9f 


:  -(2n-V 


,-1/2  lf.~l/2 


f=i  . 

pl/2  =  a.>/2  -1 


^  "  \2  ”  r  .  (c-8) 

\ 


The  meaning  and  interpretation  of  thesb  equations  are  given  in  Section  3  of  the 
main  text. 


i 


I 


I 
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The  probability  of  false  alarm  is  obtained  by  setting  R  =  0; 

Pp  =exp(-A^)  . 

t^gain,  the  signal  and  noise  spectra  do  not  enter  (C-12)  and  (C-13). 


{C-IS) 
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Appendix  D 

ASYMPTOTIC  BEHAVIOR  OF  CHARACTERISTIC  FUNCTION 
AND  ERROR  BOUND 


The  cumulative  distribution  function  (CDF)  canbe  found  from  the  character¬ 
istic  function  (CF)  by  use  of  Ref.  IS,  Eq.  (3)  : 

p(X)=l--L7.^Im{£(aexp(-ieX)l  .  (D-1) 

IT  O  ^ 

(Alternatively,  Ref,  16,  Eq.  (14)  or  (15),  could  be  used  for  the  problems  en¬ 
countered  here. )  Since  numerical  evaluation  of  this  integral  must  be  limited  to 
a  finite  interval,  it  is  necessary  to  evaluate  the  error  incurred  by  terminating 
the  range  of  the  integral.  To  evaluate  this  error,  the  behavior  of  the  CF  f(|)  for 
large  I  must  be  known. 


For  the  low-frequency  (LF)  case,  the  pertinent  CF  is  given  in  (B-31).  It 
may  be  shown  that  for  large  ?  the  second  exponential  is  upper-bounded  by  unity, 
and  the  leading  terms  behave  as 


where  u  =  ?/j3.  As  an  example  of  the  accuracy  of  this  asymptotic  form,  for 
/3  =  i  and  k  =  250,  the  exact  value  of  the  CF  is  1.489S62  •  10"’  exp  (  42.692384), 
whereas  the  asymptotic  expansion  in  (D-2)  gives  1.490021  •  10"’  exp  (-i2. 693924), 
An  upper  bound  on  the  CF  for  large  |  foilov/s  from  (D-2)  as 

(■7)  exp  [^(l-u^^)]  .  (]r)_3) 


The  error  in  terminating  the  integral  in  (D-1)  can  then  be  related  to  the  integral 

1/4 


_Lf±Lf±\  \xp[;S(i-u>/2)] 

ff  U,  n  \u/ 


>7/4 


exi 


p(^) 


r  dt  exp(-/3t)  <^— 
“'.1/2  ^ 


7/4  expEjSil-u*''^;] 


(0-4) 


T 


3/4 
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Values  of  used  to  keep  the  error  bel(w  specified  tolerances  be  found 
easily  from  (D-4).  Tliese  values  are  dependent  on  exact  value  of  d;  however, 
they  are  independent  of  Re  ^4  R?  srej  ^ere^  re,  conservative. 
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Appendix  E 

DETECTION  PROB,ABILITy  UNDER  GAUSSIAN  ASSUMPTION 

In  this  appendix^,  we  assume  the  time-bandwidth  product  large  enou^  so 
that  the  random  variable  p  in  (B-4)  can  be  asstimed  Gaussiam  In  that  event, 
we  need  only  compute  the  mean  and  variance  of  p.  We  have 


Pc  ,  l+R 
P  - —  +  1  - ; — 3-  • 

2  1+R(l_f) 

8  a 

using  (B-2),  (B-19),  (B-21),  and  (19).  *  The  variance  of  p  is  given  by 

<^2(p)  - - ff  dt  du  [4t{t)t(a)  R(t-u)  +  2R2(t-u)] 

T2  t 

.  2fR  2/3 - 1  +  exp  (-2/3)  4/3 - 1  exp  (-4j8) 

l+R(l-f)  /32  4/32 


(E-1) 


(E-2) 


using  (B-23)  and  (B-24)  for  the  case  of  interest  here.  Then  from  (B-3),  (12),  * 
and  (C-3), 


Prob(t  >X)  ---Prol.  |p  > - 2_ 

\  1+R. 


/  p  A  ' 

\<Kp)  (l+R(l-f))c(p)j 


1+R  -  A. 


[l+R(l  -f)]2/2  ri+R+3Rf - L(i +R+7Rf)  +iH  exp (-2/3)  (-4/3)1 

L  4/3  jS  4/8  J 


rhe  probability  of  false  alarm  is  obtained  by  setting  R  =  0: 


I  -  exp  (-4/3) 


/3 


^This  equation  is  in  the  main  tert  of  this  report 
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